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Abstract: We show how the Hopf algebra structure of multiple polylogarithms can be 
used to simplify complicated expressions for multi-loop amplitudes in perturbative quantum 
Q I field theory and we argue that, unlike the recently popularized symbol-based approach, 

£T) [ the coproduct incorporates information about the £ values. We illustrate our approach 

by rewriting the two-loop helicity amplitudes for a Higgs boson plus three gluons in a 
simplified and compact form involving only classical polylogarithms. 
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1. Introduction 

Higher-order quantum corrections to physical observables in perturbative quantum field 
theory require the evaluation of so-called Feynman integrals, which arise from the integra- 
tion over the momenta of (unobserved) quanta exchanged in a physical process. For this 
reason, analytical results for Feynman integrals are not only interesting in their own right 
but are also of phenomenological relevance in order to make precise predictions for current 
and future collider experiments. 



In many cases Feynman integrals can be expressed in terms of classical polylogarithms 
and Nielsen polylogarithms [1]. In the late nineties it was realized however that for certain 
multi-loop integrals new classes of transcendental functions arise that can no longer be 
expressed in terms of the classical polylogarithm functions, e.g., [2, 3, 4, 5, 6, 7, 8, 9, 10, 
11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 
35, 36]. It was soon realized that many of the new functions appearing in computations 
at higher loop orders are in fact special cases of a more general class of functions going 
under the name of multiple polylogarithms in the mathematical literature [37, 38]. While 
multiple polylogarithms do not account for all the classes of functions that can appear 
in Feynman integral computations (in some cases elliptic functions were encountered, see 
e.g., Ref. [39]), they are assumed to cover large classes of phenomenologically interesting 
Feynman integrals. 

Just like their classical analogues, not all multiple polylogarithms are independent, but 
they satisfy (complicated) relations among themselves. These functional equations make 
it possible to express a given combination of polylogarithms in a multitude of ways with 
increasing complexity. Thus, while simple and compact analytical results for Feynman 
integrals are highly desirable, the simplicity of the result might be hidden behind a swath 
of functional equations. Therefore, a systematic approach and an organizing principle 
to deal with the functional equations governing the combinatorial structure of (multiple) 
polylogarithms is a valuable tool to study Feynman integrals at higher loop orders. 

In physics a first step in this direction was made in Ref. [40] with the simplification 
of the six-point remainder function in J\f = 4 Super Yang- Mills computed in Ref. [41, 42]. 
The main tool used to achieve this simplification were the so-called symbols [58], which 
provide a way to map the combinatorics of functional equations among iterated integrals 
to the combinatorics of a certain tensor algebra. Since their introduction, symbols have 
seen many applications in physics [44, 45, 46, 47, 48, 49, 50, 51, 52, 53]. In particular, 
by now the symbols of all two-loop MHV remainder functions [54] and of the two-loop 
six-point NMHV amplitude [55] are completely known, while at three loops the symbols of 
the remainder functions for the hexagon in general kinematics [56] and for the octagon in 
special kinematics [46] are known up to some free parameters that could not be fixed from 
general considerations. However, only in the latter octagon case an integrated form of the 
symbol is also known. Indeed, while the computation of the symbol of a transcendental 
function is a straightforward and algorithmic procedure, the inverse problem (sometimes 
called integration of a symbol) of finding in an algorithmic way a function that matches 
a given symbol is currently still open, despite the fact that advances have been made in 
how to determine a class of functions that can reproduce a given symbol [57]. One of the 
problems one encounters during the integration procedure is that the symbol map is non 
injective. In particular, all terms proportional to (multiple) £ value or powers of tt are 
mapped to zero by the symbol. 

The aim of this paper is to present an alternative to, or rather an extension of, the 
naive symbol approach used in physics so far. The cornerstone of our approach is the 
coproduct on multiple polylogarithms introduced by Goncharov in Ref. [58], augmented by 
some ideas introduced in a recent paper by Brown [59]. The coproduct has the advantage 



that it does not lose as much information as the symbol, while still reproducing the latter in 
a specific limit. In particular, (multiple) ( values are not mapped to zero by the coproduct, 
thus providing valuable additional information about the function, information that could 
not be provided by the symbol alone. While we fall short of a full proof of some of our 
claims, we present accumulating evidence that our approach works in practice by applying 
it to several non-trivial cases. In particular, we consider the two-loop helicity amplitudes 
for a Higgs boson plus three gluons computed in Refs. [60, 61], where the results have been 
expressed as a complicated combination of two-dimensional harmonic polylogarithms. We 
show that the results of Refs. [60, 61] can be rewritten in a compact form that only involves 
classical polylogarithms with simple rational functions as arguments. This result confirms 
and extends a similar observation made in Ref . [53] , where it was shown that the symbol of 
the weight-four leading-color contribution of the two-loop helicity amplitudes for a Higgs 
boson plus three gluons matches the symbol of the form factor of three gluons computed 
in planar N = 4 Super Yang-Mills. 

The outline of the paper is as follows: In Sections 2 and 3 we provide a short review of 
our main actors, multiple polylogarithms and symbols. In Section 4 we give a pedestrian 
introduction to some concepts of modern algebra that are put into action in Section 5 
where we show how the Hopf algebra of Ref. [58] can be used to simplify complicated 
expressions involving multiple polylogarithms. This new technique is illustrated on some 
simple examples in Section 6, while in Section 7 we apply our tool to rewrite the helicity 
amplitudes for a Higgs boson plus three gluons computed in Refs. [60, 61] in a simplified 
form. In Section 8 we draw our conclusions. The appendices contain some technical results 
omitted throughout the main text. 

2. Multiple polylogarithms 

Analytical results for Feynman integrals are very often expressed in terms of logarithms 
and (classical) polylogarithms, 

f z dt f z dt 

lnz = / — and Li n (z) = / — Li„_i(i) , (2.1) 

Ji t Jo t 

with Lii(z) = — ln(l — z). While these functions are sufficient to describe large classes of 
Feynman integrals, it is known that especially multi-loop multi-scale integrals can give rise 
to new classes of functions. Among these new classes of functions are the so-called multiple 
polylogarithms, a multi- variable extension of Eq. (2.1) defined recursively via the iterated 

integral [37, 38] 

f z dt 

G(ai,...,a n ;z) = / G(a 2 , ■ ■ ■ ,a n ;t) , (2.2) 

Jo t - Oi 

with G{z) = 1 and where cn,z € C (they can be either constants or variables in the 

following). In the special case where all the ctj's are zero, we define, using the obvious 

vector notation a n = (a, . . . ,a), 

G(6 n ;z) = ^ln n z. (2.3) 



The number n of elements a^, counted with multiplicities, is called the weight of the multiple 
polylogarithm. 

Iterated integrals form a shuffle algebra [62], which allows one to express the product 
of two multiple polylogarithms of weight m and n^ as a linear combination with integer 
coefficients of multiple polylogarithms of weight n\ + nz, 

G(ai, . . . ,a ni ;z) G(a ni+ i, . . . ,a ni+n2 ;z) = 2_^ G(a a ^, . . . ,a ff ^ ni+n2 y,z), (2 A) 

<rG£(ni,ri2) 

where £(m,n2) denotes the set of all shuffles of n\ + rii elements, i.e., the subset of the 
symmetric group S ni+n2 defined by 

S(ni,n 2 ) = {a £ S ni+m \ <r -1 (l) < . . . < 0- -1 (ni) and <x -1 (ni + l) < . . . < cr _1 (ni+n2)} . 

(2.5) 
Whenever they converge, multiple polylogarithms can equally well be represented as 
multiple nested sums (e.g., for \xi\ < 1) [37] 

x n ^ x n 2... x n k °° x n k n k~ l "2-1 ni 

Li mi ,..., mfc (xi, . . . ,x k ) = 2_^ ZjmZJm ZW = / > ZW / J ■■■ / J ZW • 

ni<ri2<--<n k i z k nfc=l K 7ifc_i=l m=l - 1 - 

(2.6) 
Note that we are using Goncharov's original summation convention [37]; other authors 
define Li mii ... ;mfc (xi, . . . , x^) using the reverse summation convention instead, i.e. ni > 
• • • > nfc. The G and Li functions define in fact the same class of functions and are related 
by 

Li mi ,..., mfe (xi, . . . ,x k ) = (-l) k G( 0, . . . , 0, -1 . . . , 0, . . . , 0, 1 - ) . (2.7) 

jrtj, — 1 mi — 1 

It is possible to find closed expressions for special classes of multiple polylogarithms 
in terms of classical polylogarithm functions, e.g., for o/Owe have 

G(0 n ; z) = -. In™ z, G(a„; z) = - In" (l - -) , 
G(0 n _i,a;z) = -Li n ^-J , G(0„,a p ;z) = (-1) P 5„ )P ^-J , 

where S^p denotes the Nielsen polylogarithm [1]. 

All the notations so far follow the conventions used in physics. Some of the formulas 
used later on however take a nicer form in a different notation commonly used in the 
mathematical literature, 

ra n+1 fo 

I(a ;ai,. . . ,a n ;a n+ i) = — — I(a ;ai, . . . ,a n -i;t) . (2.9) 

Jao t a n 

The notations (2.2) and (2.9) are related by (note the reversal of the arguments) 

G(a n ,... ,ai;a n+ i) = 1(0; ai,. . . ,a n ;a n+ i) . ( 2 -10) 

The iterated integrals defined in Eq. (2.9) are slightly more general than the ones usually 
defined by physicists, as they allow to freely choose the base point of the integration. It is 



nevertheless easy to convert every integral with a generic base point ao into a combination 
of iterated integrals with base point 0. An example will clarify this. First, it is easy to see 
that at weight one we have 

I(a ; ai; a 2 ) = 1(0; au a 2 ) - 1(0; oi; a ) = G(oi; o 2 ) - G(oi; a ) • (2.11) 

Starting from weight two the relation is more complicated because of the nestedness of the 
integration, e.g., 

f as dt f a ' A dt 

I(ao;ai,a 2 ;a 3 ) = - I(a ;ai;t) = [l(0;ai;t) - J(0;ai;a )] 

Jao t-a 2 J ao t-a 2 

(2 12) 
= J(0;ai,a 2 ;a 3 ) - l(0;ai,a 2 ;a ) - l(0;ai;a )[i"(0;a 2 ;a 3 ) - l(0;a 2 ;a )} 

= G(a 2 ,ai;a 3 ) - G(o 2 ,oi; a ) - G(oi; a )[G(a 2 ; a 3 ) - G(a 2 ;a )] . 

In the rest of the paper we mostly use the '/' notation used in the mathematical literature, 
as it makes most of the formulas much simpler, keeping in mind that we can easily recover 
the 'G" notation via the above procedure. 

Just like their classical analogues, multiple polylogarithms satisfy a large number of 
functional equations among themselves. When expressing a Feynman integral in terms of 
multiple polylogarithms, we can therefore arrive at a complicated combination of multi- 
ple polylogarithms, which, if the corresponding functional equations were known, could 
potentially be reduced to a much simpler expression. While these functional equations 
are unknown in general, they can often be circumvented in practice by using the so-called 
symbol, which we will review in the next section. 

3. Symbols 

In this section we give a short review of symbols [43]. Symbols were first introduced in 
physics in Ref. [40] where they were used to simplify the six-point remainder function in 
J\f = 4 Super Yang-Mills computed in Ref. [41, 42]. The main idea is to map a (compli- 
cated) combination of multiple polylogarithms to a certain tensor algebra over the group 
of rational functions (the tensors being called symbols) such that, at least conjecturally, all 
the functional equations among the polylogarithms are mapped to simple algebraic iden- 
tities in the tensor algebra. Currently, two different definitions of symbols are in use in 
physics, 

1. In Ref. [40] the symbol of a transcendental function F w (x±, . . . ,x n ) of weight w in 
the variables x±,. . . ,x n was defined recursively by considering the total differential 
of the function F w . More precisely, if the total differential of F w can be written in 



the form 



dF w = J2 F i,w-idlnRi, (3.1) 



where Fi W —i are transcendental functions of weight w — 1 and the Ri are rational 
functions in the variables x±, . . . , x n , then the symbol of F w can be computed recur- 
sively in the weight by 

S(F w ) = YS(F iyW _ 1 )®R i . (3.2) 






Multiple poly logarithms indeed satisfy a differential equation of the type (3.1) [63], 

dl(a ;ai,... ,a n ;a n+1 ) = y~]l(a ;ai,... ,6*,... , a n ; a n+1 ) din I — ) , (3.3) 

i=i ' 

where the hat indicates that the corresponding element is omitted. We emphasize 
though that Eq. (3.3) is strictly speaking only valid if all the at are generic, i.e., non 
zero and mutually different. In the non-generic case the differential equation (3.3) 
can take a different form [63]. 

2. In Ref. [57] an alternative definition of a symbol was given. It was shown that the 
symbol of a multiple polylogarithm can be obtained by summing over certain dissec- 
tions of a rooted and decorated polygon associated to a multiple polylogarithm [64], 
and the combinatorics of these dissections reproduces precisely the symbol obtained 
from the recursive procedure of Ref. [40] . 

Both definitions have their advantages and disadvantages. While the recursive definition 
has the obvious advantage that it is not necessarily restricted to multiple polylogarithms 
but extends to any class of (transcendental) functions defined through iterated integrals and 
satisfying a differential equation of the type (3.1), the second definition maps the combina- 
torics of the symbol to the combinatorics of rooted decorated polygons, a correspondence 
currently only established in the case of polylogarithms. On the other hand, the approach 
based on polygons is algebraic in nature, and does not make any difference between con- 
stants and variables. As an example, the differential equation approach would assign a zero 
symbol to In 2 (as <iln2 = 0) while <S(ln2) = 2 from the polygon approach. As otherwise 
both definitions are equivalent and give the same answer for multiple polylogarithms, we 
will in the following not distinguish them any further. 

The symbol map S fulfills various properties. First, S is linear and maps a product 
of multiple polylogarithms to the shuffle product of their tensors (more precisely, S is an 
algebra homomorphism, see Section 4). Next, each factor in the symbol is additive with 
respect to multiplication, 

...<g)(a-6)(8>... = ...<8)a<g>... + ...<8>&<8>.... (3.4) 

This implies in particular that 

...01®... = 0, (3.5) 

and more generally if p n denotes an n-th root of unity, 

...<g>p n <g>... = 0. (3.6) 

From Eq. (3.4) and Eq. (3.5) it is clear that each factor in the tensor product 'behaves as 
a logarithm'. 

The first and the last entry of the symbol of a function carry some special information. 
Let us consider a transcendental function F w (xi, . . . , x n ) whose symbol takes the form 

S(F w (xi,...,x n )) = ^2 c ii,...,i w Uh®---®Ui w , (3.7) 



where Cj 1 j^ are (rational) numbers and uji k are rational functions in the X{. The symbol 
of the derivative of F w is given by 

S[—F w {xi,...,x n )\= ^2 c hj ... j i w u>i 1 ®...®uji w _ 1 (—lnuji w j. (3.8) 

In other words, a derivative only acts on the last entry of the symbol. 

The first entry of a symbol encodes in a similar way the information about the mon- 
odromies (discontinuities) of the function F w . More precisely, if A4 Xk=a is the operator 
that computes the monodromy of F w around x^ = a, then 

S(M Xk=a F w (xi,...,x n )) = ^2 (M^alnu;^) Ci u ... tiw w i2 ® ...®uj iw . (3.9) 

Note that the action of the monodromy operator is trivial on the left-hand side, because 
it only acts on ordinary logarithms, 

, . , I 2-7ri , if Ui-, has a zero for x^ = a , ,„ _, n . 

M Xk= alnuJ h = < n ,, . (3.10) 

I , otherwise . 

We prefer nevertheless to write Eq. (3.9) in this apparently more complicated form in order 
to exhibit the duality to Eq. (3.8). 

So far we have only dealt with the problem of how to compute the symbol of a func- 
tion. Indeed, using any of the two definitions we can compute the symbol of any linear 
combination of products of multiple polylogarithms. Once the symbol has been obtained, 
the identities (3.4) and (3.5) allow us to simplify the symbol, which is equivalent to ap- 
plying functional equations to the original expression. We then have to face the problem, 
however, of finding a (simpler) function with the same symbol. While there are rules how 
to compute the symbol of any given combination of polylogarithms, the inverse step of 
integrating the symbol to a function (i.e., of finding a function with the same symbol) is in 
general much more difficult. In Ref. [57] a prescription was given that allows one to make 
an educated guess for the class of functions that can give rise to a given symbol. Once 
such a class of functions has been determined, one can write down a linear combination 
(with some free coefficients) of these functions and equate their symbols, obtaining in this 
way a linear system for the coefficients. However, even after this step there is a remaining 
ambiguity because the symbol map is not injective. As an example, we have 

S(iir)=0 and S(( n ) = 0. (3.11) 

As S maps products of functions to shuffle products of tensors, Eq. (3.11) implies that all 
terms proportional to ( values and / or i?r will be mapped to zero by S. As a consequence, 
even if we succeed in finding a simpler function with the same symbol as our original 
function, we are unable to fix the terms proportional to, e.g., £ values based on the symbol 
alone. 

The aim of this paper is to introduce a framework similar in spirit to the symbol, but 
where terms proportional to (, values and m are not mapped to zero. Such a framework 



should of course retain all the good features of the symbol and reduce to the definition 
of a symbol in a suitable limit. In the rest of this paper we argue that such a framework 
is provided by the Hopf algebra of multiple polylogarithms introduced by Goncharov in 
Ref. [58], augmented by some ideas inspired by a recent paper by Brown [59]. 

4. Algebras, coalgebras and Hopf algebras 

The aim of this section is to provide a (pedagogical) review of the algebraic notions used 
throughout the rest of the paper. The content of this section is standard textbook material 
in mathematics. We nevertheless include it here because, at least to our knowledge, most of 
these concepts have only been rarely used in the context of Feynman integral computations. 
We emphasize that we do not aim at providing a rigorous mathematical exposition of the 
different topics, but rather content ourselves to provide a pedestrian introduction, and 
we refer to the dedicated mathematical literature for further details. In particular, we 
will proceed by analogy with similar mathematical concepts that are of everyday use in 
physics. Furthermore, we will not be concerned about technical details, such as for example 
the differences between rings and modules on the one hand, and vector spaces and fields 
on the other hand. As a consequence, we will use the different notions interchangeably in 
the following. 

4.1 Algebras: first definition 

We start by reviewing some basic notions about algebras. An algebra over a field K (= R 
or C in general) is a K- vector space A together with a map 



m : Ax A ->■ A 

(a, b) i->- m(a, b) = a • b 



(4.1) 



that is associative 



and has a unit e, 



c) = (a-b)-c, (4.2) 



e ■ a = a • e = a . (4-3) 

Furthermore, the multiplication is compatible with the vector space structure, i.e., dis- 
tributive, 

a-{b + c) = a-b + a-c and (a + b)-c = a-c + b-c, (4-4) 

and associative with respect to scalars, 

a-(kb) = k(a-b) and k (£a) = (k£) ■ a . (4.5) 

where k,£ € K are scalars. Let us highlight at this stage some features that will be 
useful later on. First, the distributivity (4.4) implies that the multiplication m is in fact 
a bilinear map from A x A to A. Second, as a consequence of the compatibility with the 
scalar multiplication, Eq. (4.5), we will assume from now on that the field of scalars K is 



part of the algebra itself, i.e., IK can be embedded into A. Under this assumption it is easy 
to see that we can identify the unit element e in Eq. (4.3) with the unit element 1 G IK. 

While this definition of an algebra is presumably familiar to most physicists, it is not 
well suited to understand the link to coalgebras and Hopf algebras. For this reason, we 
will now reformulate the above definition in terms of tensor products of vector spaces. 

4.2 Tensor products of vector spaces 

As the concept of tensor product used in the following is different from the definition 
commonly used in physics, let us review the mathematical definition of a tensor product. 
Consider three vector spaces U, V and W. A standard textbook result then states 
that there is a unique vector space, called the tensor product U <8) V of U and V, and a 
unique bilinear map t:UxV—>U®V such that for every bilinear map ft : U x V — >■ W 
there is a unique linear map fx such that 

P = HT, (4.6) 

where /it denotes the composition of jji and r. The bilinear map r simply assigns to two 
vectors a and b their tensor product, i.e., r(a, b) = a®b. In other words, we can reformulate 
Eq. (4.6) in an equivalent, but more accessible, form: for every bilinear map j3 there is a 
unique linear map [i such that 

/3(o, b) = fi(a <g> b) . (4.7) 

According to Eq. (4.7) it is always possible to break the bilinearity of /3 up into the bilin- 
earity of the tensor product and the linearity of fi. As an example, we have 

P(k a + £ b, c) = fi((k a + £ b) <g> c) = fi(k (a <g> c) + £ (b <g> c)) 

(4.8) 
= k fi(a ® c) + £ fi(b ® c) = k /3(a, c) + £ /3(b, c) . 

Before applying this result to the definition of algebras, let us take the opportunity to 
introduce a diagrammatic tool useful to describe algebraic structures. If we represent each 
map between vector spaces by an arrow, e.g., 

U xV ^U®V, UxV^W, U®V^W, (4.9) 

then the relation (4.6) between r, j3 and fj, is conveniently described by the following 
commutative diagram, 

U xV U®V 



/? 



H (4.10) 



W 

The word 'commutative' refers in this context to the fact that we can take any path through 
the diagram from U x V to W and we always arrive at the same result. This is precisely 
the content of Eq. (4.6). 



4.3 Algebras: second definition 

We have seen that the distributivity of the multiplication m implies that m is in fact a 
bilinear map from A X A to A- Thus, following our considerations on tensor products, 
there is a unique linear map fj, : A ® A — )■ A such that 



a • b = m(a, b) = /j,(a ® b) . 



(4.11) 



The associativity of the tensor product can then be summarized by the following property 
of the linear map //, 

//(id <g> n) = n(n ® id) . (4.12) 

Indeed, acting with ^(id ® /i) on a ® b ® c gives 1 

/Lt(id <g> /i)(a <g 6 <g> c) = //(a (8 //(& (S) c)) = /i(a (8) (6 • c)) = a • (6 • c) , (4.13) 

while acting with /i(/i® id) gives 

(jl(h id) (a (g) 6 (g) c) = /u(/i(a (g 6) (g> c) = //((a ■ 6) <g c) = (a ■ 6) ■ c. (4.14) 

In order to make the transition to coalgebras easier, it is useful to rephrase Eq. (4.12) in 
terms of commutative diagrams. It is easy to check that Eq. (4.12) is equivalent to the 
commutativity of the diagram 



A® A® A 



id <g n 



fj, (g id 



A® A 



ft 



A® A 



/' 



(4.15) 



A 



The existence of a unit element can also be recast into this new language. Indeed, we 
previously assumed that the field of scalars is embedded into the algebra. This implies 
the existence of a map e : K — > A that assigns to each scalar k and element e(k) £ A. 
Compatibility with the scalar multiplication forces us to require that e be linear and that, 
for any vector a, multiplication with the scalar k or the vector e{k) gives the same result, 



e(k)-a = ka and e(k) ■ e(£) = e(k £) . 

In terms of commutative diagrams Eq. (4.16) is equivalently expressed as 

K®A 



e<gid 



(4.16) 



A® A 




(4.17) 



For two functions / and g we define (/ ® g) (a ® b) — f(d) (S> 9(b)- 
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where s denotes the scalar multiplication of a vector and a scalar, s(k, a) = ka. Note that 
the embedding e, together with Eq. (4.16), implies the existence of the unit e = e(l) in the 
usual sense, because 

e(l)-o = lo = o. (4.18) 

For this reason, the map e is usually referred to as the unit of the algebra. 

Another recurrent theme in the study of algebraic structures is the study of the 
structure-preserving maps (the so-called homomorphisms) . If A and B are algebras, then 
a homomorphism from A to B is a linear map (f> that preserves the multiplication, 

4>{a ■ b) = <j)(a) ■ (f>{b) . (4.19) 

At this stage we can identify the multiple polylogarithms as an algebra: the multipli- 
cation is given by the shuffle product, while the scalars are given by (rational) numbers. 
In addition, the shuffle product preserves the weight (the product of two polylogarithms 
of weight ri\ and n<i gives a linear combination of polylogarithms of weight n± + ri2). This 
feature is formalized by the notion of a graded algebra, i.e., an algebra that is a direct sum 
as a vector space, 

oo 

-4 = A*, (4.20) 

?1=0 

such that the multiplication preserves the grading, 

An element of A n is said to be of weight or grade n. The multiple polylogarithm algebra is 
thus graded by the weight. In the following we always assume that the weight part of a 
graded algebra coincides with the field of scalars, Aq = K. Note that this is in agreement 
with our naive expectation in the case of multiple polylogarithms: the weight part of the 
algebra consists of all objects of transcendental weight 0, e.g., rational numbers or rational 
functions. 

To summarize, we arrive at the following conclusion: an algebra A is a vector space 
equipped with a linear map \x : A C3 A — > A and an algebra homomorphism e : IK — >■ A such 
that /i satisfies the associativity condition expressed through Eq. (4.12) or, equivalently, 
through the commutative diagram (4.15). The multiple polylogarithms, equipped with the 
shuffle product, then form an algebra graded by the weight. 

So far we have only formalized the algebraic structure of multiple polylogarithms in an 
unusual way. In Section 5 we will see that multiple polylogarithms carry another structure, 
a so-called Hopf algebra structure, which is the topic of the rest of this section. 

4.4 Coalgebras as the duals of algebras 

In this section we introduce coalgebras and coproducts, one of the main topics in the rest 
of this paper. We refrain from giving a mathematically precise axiomatic definition of 
a coalgebra, and we rather proceed by analogy with similar concepts of everyday use in 
physics: vector spaces and hermitian conjugates of linear maps. 
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Let us start by reviewing some textbook material on the duals of vector spaces. Con- 
sider two vector spaces V and W and denote their duals 2 , i.e., the vector spaces of all 
linear forms on V and W, by V* and W*. Let A be a linear map from V to W. It is 
well-known that A induces a linear map between the duals, the hermitian conjugate A\ in 
the opposite direction, 

A. TTr , Tlr * A\ 



V 



W and W* 



V* . 



(4.22) 



Eq. (4.22) provides us with a diagrammatic rule to derive the commutative diagram describ- 
ing some algebraic structure when replacing all the vector spaces by their duals. Indeed, 
when dualizing, we 

1. replace each vector space by its dual, 

2. replace each linear map by its dual (= hermitian conjugate), 

3. reverse all the arrows in the diagram. 

If we upgrade V to an algebra, then it comes equipped with two linear maps, the multipli- 
cation jjl and the unit e, so we can consider their 'hermitian conjugates', i.e., the structures 
induced on the dual space V* by the multiplication and the unit. 

A coalgebra is defined as the dual of an algebra, i.e., if A is an algebra with multipli- 
cation jjl : A <S> A — >■ A and unit e : IK — )■ A, its dual C = A* is equipped with linear maps 
A = /v : C — > C®C (the comultiplication) and rj = e' : C — > K. (the counit). The properties 
of the algebra operations (associativity and unit) will also be reflected in the coalgebra 
C. Using the diagrammatic rule given in Eq. (4.22), we can easily obtain the commuta- 
tive diagrams that describe the duals to the associativity (4.15) and the unit (4.17) (the 
coassociativity and the counit), 



c®c®c 



A <8) id 



id® A 



C®C 

A 



C®C 



A 



r\ (g> id 




(4.23) 



The coassociativity can also be written in equations as 

(id <g> A)A = (A <g> id) A . 



(4.24) 



As the coassociativity will be extensively used in the remainder of this paper, let us elabo- 
rate on it for a while. While the multiplication fi corresponds to the operation of 'multiply- 
ing together' and the associativity expresses the fact that it is immaterial in which order 
we multiply three or more elements together, the comultiplication A morally corresponds 
to the operation of 'decomposing', and the coassociativity asserts that the order in which 



2 In quantum mechanics, the elements of V are called 'ket's, while the elements of V* are called 'bra's 
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this operation is iterated is irrelevant. To be more concrete, take an element a G C, and 
consider its coproduct, schematically written as 

A(a) = ^af ® af } , (4.25) 

i 

for some af € C. By acting with A on a, we have decomposed it into (a combination 
of) two pieces. We can iterate this process, and decompose Eq. (4.25) further into (a 
combination of ) three pieces. At this stage, we have the choice to decompose either a\ 

(2) 

or a\ , 

A K (1) ) = E <$f ® <? and *(«i 2) ) = E «g 1} ® «if ( 4 - 26 ) 

In the first case we arrive at 



(a ® id)A(o) = Yl 4T ® «S 2) ® «f ) > ( 4 - 27 ) 



(id ® A)A(a) = J^ af ® ag X) ® ag 2) . (4.28) 



^C 1 . 1 ) ^ „(L2) ^ „(2) 

while in the second case we arrive at 

,(1) ,,1.2.1) J 2.2! 
i 



While Eqs. (4.27) and (4.28) are in general different, coassociativity asserts that the two 
expressions must be equal 



(a ® id) a(o) = J2 a ff ® °S 2) ® a *- 2) = E a i 1] ® °2 ;1) ® a S' 2) = ( id ® A ) A ( a ) • ( 4 - 29 ) 



As a consequence, there is an essentially unique way to iterate the coproduct, 

H A^®?^ d ^®^®^ A ^f d ... (4.30) 

4.5 Bialgebras and Hopf algebras 

We are now only one step away from defining a Hopf algebra. First, a bialgebra is an 
algebra that is at the same time a coalgebra, i.e., a vector space equipped with both a 
multiplication fx and a comultiplication A 3 . We emphasize that in this setting [i and A 
are independent and in general not hermitian conjugate to each other. Furthermore, we 
require the multiplication and the comultiplication to be compatible with each other, i.e., 
the coproduct of a product equals the product of the coproducts (in other words, A is an 
algebra homomorphism) , 

A(o • b) = A(o) • A(6) , (4.31) 

where the multiplication in the right-hand side is taken in each factor of the tensor product 
separately, 

(ai ® a 2 ) • (6i ® b 2 ) = (oi • &i) ® (a 2 • b 2 ) ■ (4.32) 



3 We will from now always tacitly assume that all (co)multiplications are (co)associative. 



13 



A Hopf algebra % is a bialgebra equipped with an additional structure, the so-called 
antipode S :H — > rl satisfying the properties 

S(a ■ b) = 5(6) ■ S(a) and //(id ® 5)A = n(S ® id) A = . (4.33) 

As in the rest of this paper we do not make explicit use of the antipode, we do not elaborate 
on it any further. 

Let us conclude this section by introducing some notations that will be useful in sub- 
sequent sections. Consider a Hopf algebra % with coproduct A, and assume that % is 
graded (as will be the case for the multiple polylogarithms) , 



H = @'H n . (4.34) 



n=0 

If the coproduct respects the weight, we can decompose the action of the coproduct ac- 
cording to 

n n A H p ®H q . (4.35) 

p+q=n 

We can then write the action of A on T~i n as 

A= Y, A ^ ( 4 -36) 

p+q=n 

where A p ^ q is the part of the coproduct that takes values in % p <g> H q . In a similar way we 
define A pA ... ]T . as the component of the iterated coproduct that takes values in H p ® H q <g> 
. . . <8> Tir- Finally, it is sometimes useful to define the reduced coproduct A' via 

A(o) = 1<8)o + o<8)1 + A'(a) . (4.37) 

An element a G % such that A' (a) = is called a primitive element of %. 

5. The multiple polylogarithm Hopf algebra 

In this section we apply the algebraic concepts of the previous section to multiple polylog- 
arithms. As a result, we obtain a framework that contains the symbol in a certain limit, 
but is more general and incorporates, in particular, the Q values. 

As a starting point, let us denote by % the algebra formed by the multiple polylog- 
arithms equipped with the shuffle product. We already know that H is graded by the 
weight of the polylogarithms. In Ref. [58] Goncharov showed that T-L can be equipped with 
a coproduct which turns it into a Hopf algebra. The coproduct on multiple polylogarithms 
is given by [58] 

A(J(ao;ai,... ,a„;a„+i)) 

" k 

Y[I{ai p ;a ip+ i,... ,a ip+1 -i;a ip+1 ) 

.p=0 

(5.1) 



/ J I(ao',a>h, ■ ■ ■ >Oj fc ;a n+ i) 

0=Ji<J2<---<*fe<*fc+i=n. 
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The fact that Eq. (5.1) defines a genuine coproduct, i.e., that A is coassociative, Eq. (4.24), 
and an algebra homomorphism, Eq. (4.31), is a non-trivial statement. In addition, Eq. (5.1) 
preserves the weight, i.e., the sum of the weights in each term is equal to n. We stress 
that Eq. (5.1) is strictly speaking only valid when all the ctj's are generic, i.e., non zero and 
mutually different. The definition of the coproduct in the non- generic case involves several 
technical steps that do not add anything new to the discussion in the main text of the 
paper, and we refer to Appendix A or to Refs. [58, 63] for the definition of the coproduct 
in the non-generic case. Let us quote here only the explicit formulas for the coproducts for 
the ordinary logarithm and the classical polylogarithm, 

A(ln z) =l<g)lnz + lnz<g>l, 

n ~ 1 l k (5 2) 

A(Li n (z)) = 1 <g> U n (z) + Li n (z) ® 1 + J2 u n-k(z) ® —r£- . 

fe=i 

Eq. (5.2) is enough to compute the coproduct of any expression made out of ordinary 
logarithms and classical polylogarithms only. Indeed, we can use Eq. (4.31) to obtain for 
example, 



A (In x In y) = A (In x) A(ln y) = [1 ® In x + In x ® 1] [1 ® In y + In y ® 1] 
= 1 ® (In x In y) + In x ® In y + In y ® In x + (In x In y) ® 1 . 

Furthermore, it is easy to prove the following result, 



(5.3) 



n 

' n 



A(ln n z) = ^[ \\n k z® ln n - fe z . (5.4) 

fc=o ^ ' 

The coproduct can be used to simplify expressions involving polylogarithms in the 
same way as the symbol. Indeed, suppose that we have two expression F w and G w of 
weight w that are equal (modulo functional equations). Then it is clear that also their 
coproducts must be equal, 

A(F W ) = A(G W ) , (5.5) 

and also 

A'(F W ) = A'(G W ) . (5.6) 

It is important to note that Eq. (5.6) only involves polylogarithms of weight w' < w. As 
a consequence, it is enough to know the functional equations of lower weight in order to 
check the equality. These functional equations of lower weight might themselves still be 
complicated or unknown, so we have apparently not gained anything. In such a scenario 
we can iterate the procedure by applying the coproduct again to one of the factors in 
the tensor product, and the coassociativity of the coproduct ensures that this iteration is 
unique. In this way we obtain a whole tower of expressions, which at each stage involve 
only transcendental functions of lower weight, 

F w = G w ->• A(F W ) = A(G W ) ->• (id ® A)A(F W ) = (id ® A)A(G W ) ->• . . . (5.7) 

As an example, in the case of a function of weight four, we obtain the following identities, 
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F 4 = G 4 



A 3 ,l(*4 



A 2 ,l,l(^4) 




Ai 3 (G 4 ) 



Ai,i,2(G 4 ) 



In the extreme case where we go down to Ai i, we have decomposed a weight w poly- 
logarithm into a tensor of rank w made out of polylogarithms of weight one, i.e., ordinary 
logarithms, for which all the functional equations are known. It can be shown that in this 
way we obtain precisely the symbol of the function (up to some technical detail that will 
be discussed later). In other words, the symbol is nothing but the maximal iteration of the 
coproduct. Besides providing a precise definition of the symbol, this approach also shows 
why the symbol alone is insufficient to determine the function completely. Indeed, requiring 
two expressions to have the same symbol is equivalent to require that they give the same 
result when acted upon with Ai i. While this approach has the obvious advantage that it 
reduces the problem to the sole application of functional equations for ordinary logarithms, 
it does in general not imply the equality of the other components of the coproduct. The 
information on the terms that are missed by the symbol is nevertheless contained in these 
other components (at least to some extent). To see how this works, and to see how the 
£ values that were missed by the symbol arise in the other components, we first need to 
overcome some technical obstacles that we will discuss now. 

5.1 The multiple £ value Hopf algebra 

In the previous section we argued that the coproduct provides a more general 'calculus' 
that contains the symbol in some limit and that cures most of the unwanted features of 
the symbol. However, a consistent 'extended symbol calculus' should be compatible with 
specializations of the arguments. 

Multiple C values are, by definition, the values of the multiple polylogarithms (in the 
series representation) with all arguments equal to unity. It is easy to check that they form 
in fact a sub-Hopf algebra Z of the multiple polylogarithm Hopf algebra H. From Eq. (5.2) 
one immediately sees that the coproduct of Q values of depth one is given by, 



A(Cn) = A(Li n (l)) = 1 ® Cn + Cn ® 1 



(5.8) 
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i.e., £ value of depth one are primitive elements in Z, and thus in %. 

At this point we have to face a subtle problem for the even £ values. We know that 
the even Q values are not independent, but they are all proportional to powers of (2, e.g., 



C4 = ^C|- (5.9) 



Thus, 



A(C 4 ) = ^ A (C 2 ) 2 = ^[1 ® C2 + C2 ® I] 2 = ^[1 ® (I + Cl ® 1 + 2C 2 ® C2] , (5-10) 

and so there is a contradiction with Eq. (5.8), unless 'C2 = 0', i.e., unless we work modulo 

C2, 

A(C 2 ) = 0. (5.11) 

As a consequence, we lose all information on the terms proportional to it 2 in the coproduct. 
Hence, if this was the case we would not have gained anything over the naive symbol 
approach. 

In Ref. [59] Brown argues that instead of defining the coproduct of C2 to be zero, it is 
consistent to define 

A(C 2 ) = C2®1, (5-12) 

and more generally 

A(C2n) = C2n ® 1 • (5.13) 

This definition obviously solves the problem we had before, because 

A(C 4 ) = ^A(C 2 ) 2 = ic 2 ® I] 2 = ^ C 2 2 ® 1 = C4 ® 1 • (5-14) 

15 15 15 

Even though Eq. (5.12) was introduced in Ref. [59] in the context of multiple £ values, 
we argue that it equally well holds in more general situations. Moreover, we conjecture 
that Eq. (5.12) can be extended to 

A(vr)=7r®l. (5.15) 

This definition is obviously consistent with Eq. (5.12). In addition, it allows to extend the 
coproduct to include the m terms in a consistent way. A word of caution is however in 
order: due to the monodromy of the logarithm, we should define 2i7r®x = 0, Vx, and thus 
also 4-7T 2 ® x = 0. In practice though, we observed that we never have to worry about the 
monodromy of the logarithm in physical applications. Indeed, in a physical computation 
the Riemann sheets of the logarithms are fixed, e.g., by assigning a small imaginary part 
to x, such that 

ln(x + iSs) = In \x\ + 5m9{— x) . (5.16) 

The definition (5.15) changes the coproduct from being a map 

/\:n^%®U, (5.17) 
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to a map 

A:U^U%U* , (5.18) 

where W* is the quotient of % by the (two-sided) ideal 4 generated by it. The iterated 
coproduct then takes the form 

H A«®7r^«®«' r ®?r A ^? id .... (5.19) 

Loosely speaking, this means that we drop all powers of n in all factors of the coproduct, 
except for the first one. 

We are now able to state our main conjecture. We conjecture that using the coproduct 
together with the definition (5.15) we obtain an extension of the symbol calculus that takes 
into account the £ values as well. More precisely, if we have a function F w of weight w and 
if we can find a (simpler) function G w such that 

A'(F W ) = A'(G W ) , (5.20) 

then 

F w = G w + 2 c i Pw,i i (5-21) 

i 

where the sum runs over all primitive elements of weight w of 7i for some (rational) coeffi- 
cients Ci. We see from Eq. (5.21) that the reduced coproduct does still not entirely fix the 
function. A similar observation was made in Ref. [59] in the case of multiple £ values. In 
practice, the primitive elements turn out to be constants of a given weight, e.g., 

1. powers of it, 

2. ( values of depth one, ( n , 

3. Clausen values at the roots of unity, 

Cl n (^)=K n U n (e^ N ) , (5.22) 



N J 
where 9K n denotes the real part for n even and the imaginary part for n odd. 

Even though the function is not entirely fixed, we believe that this approach constitutes 
an important improvement over the pure symbol approach. Indeed, while a pure symbol 
approach misses for example all functions multiplied by multiple £ values, only very few 
constants are left undetermined by the coproduct. The undetermined constants can easily 
be fixed by, e.g., comparing to numerical values at a few points, requiring the function to 
have the right limits, etc. While we fall short of a full proof of our conjecture, we have 
checked the consistency of our 'extended symbol calculus' by applying it to hundreds of 
functional equations among multiple polylogarithms. A selection of results will be shown 
as an illustration in Sections 6 and 7. 



4 We recall that an ideal I in a commutative algebra A is an additive subgroup such that Va £ A and 
VV G X, we have a ■ r £X. 
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Let us conclude this section by discussing how differential and monodromy operators 
act on the coproduct, i.e., how to generalize the relations (3.8) and (3.9) to our framework. 
We conjecture that 

A (^-)-( M *^) A(F - ) ' (5.23) 

A (M Xk=a F w ) = (M Xk=a id) A(F W ) , 

i.e., differential operators only act in the last component of the coproduct, while mon- 
odromy operators only act in the first component. Note that the same statement is true 
for the iterated coproduct. While we fall short of a full proof of Eq. (5.23), we were able 
to check our claim in the special case where F w is a multiple polylogarithm with generic 
arguments. The proofs of these special cases are presented in Appendix B. 

5.2 Relationship between the coproduct and the symbol 

In this section we briefly discuss the relationship between the coproduct and the symbol. 
While it is possible to proof in general that the combinatorics of the maximal iteration 
of the coproduct on multiple polylogarithms matches precisely the combinatorics of the 
maximal dissections of the rooted and decorated polygon associated to a polylogarithm [65], 
we do not present a firm proof in this paper, but merely state the observation that this 
correspondence holds in the all cases we have considered. We only motivate the relationship 
by analyzing how the coproduct behaves under differentiation. If F w is a transcendental 
function of weight w, then without loss of generality we can write its iterated coproduct in 
the form 

Ai ) ... > i(F w ) = ^Ai ) ... ) i(F i)W _i)®liii2 i . (5.24) 

i 

If we now act with (id (g> d) on this expression, we obtain, using Eq. (5.23), 

Ai,...,i(dF w ) = J2Ai,...,i(Fi, w -i)®d\nRi, (5.25) 

i 

i.e., we obtain an expression which is dual to the differential equation (3.1) defining the 
symbol. We emphasize that we claim in no way that this provides a proof of the fact that 
the maximal iteration contains the symbol, but we hope that it provides a feeling to the 
reader why this relationship is true. Note however that the symbol is not exactly equal to 
Ai i. Indeed, the symbol does not contain any information about terms proportional to 
in, whereas these terms are incorporated into Ai i through Eq. (5.15). In other words, 
the correct relationship between the symbol and the maximal iteration of the coproduct 
reads 

S = Ai,...,! mod 7T . (5.26) 

Even though we did not provide a proof of Eq. (5.26) , we can show that it is consistent 
with our general knowledge about symbols: 

1. the fact that each entry in a symbol is additive like a logarithm, Eq. (3.4), is a 
consequence of the fact that each factor in the tensor product on the left-hand side 
of Eq. (5.26) is a logarithm. 
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2. the derivative identity in Eq. (5.23) reduces to Eq. (3.8) if we restrict ourselves to 
the maximal iteration of the coproduct. This is obvious from Eq. (5.25). 

3. similarly, the monodromy identity in Eq. (5.23) reduces to the corresponding identity 
for the symbol, Eq. (3.9). 



6. Examples 

In this section we present some simple examples of how the coproduct can be used to 
simplify expressions involving multiple polylogarithms. The examples in this section do 
not provide any new results, but they are simple enough so that all the steps can be carried 
out by hand. They are therefore rather meant to illustrate how to use the coproduct in 
practise to perform computations. 

6.1 Inversion relations 

We start by considering inversion relations for classical polylogarithms. Throughout this 
section we assume that x is a real positive variable to which we assign a small positive 
imaginary part. 

We proceed in a bootstrap and build up the inversion relations by a recursion in the 
weight. For the classical polylogarithm of weight 1, the inversion relation is easy to obtain, 



Li i I- 



ln I 1-- 

x 



ln(l — x) + ln(— x) 



ln(l — x) + lnx — in . (6.1) 



In order to obtain the inversion relation for weight 2, we act with Ai i on Li2(l/x) and 
insert the inversion relation for Lii(l/x), 



A 



1,1 



Lio 



In 1 



1 



In 



x J \x / 

ln(l — x) <S> In x — In x d$> In x + in CED In x 



(6.2) 



A 



1,1 



Li2(x) In x + in\n.x 



Following our conjecture, we conclude that the arguments on the left and right-hand sides 
are equal modulo primitive elements of weight two. We thus make the ansatz, 

1\ .. , , 1 



Li 2 



-Li 2 (ar) 



In x + in In x + c n , 



(6.3) 



for some rational number c. Specializing to x = 1, we immediately obtain c = 1/3, which 
is indeed the correct inversion relation for Li 2 . We emphasize at this stage the importance 
of the definition (5.15). 

Moving on to weight 3, we act with Ai 1 1 on Li3(l/x) and obtain 



A L 



1,1 



Li, [ - 



In ( 1-- 

x 



In ( - 
x 



In ( - 
x 



ln(l — x) <8> In x <8> In x + In x <8> In x <8> In x — in <8> In x <g> In x (6-4) 



Ai,i, 



Li 3 (x) + -ln 3 x 
o 



in 



In x 
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Eq. (6.4) is not yet the correct inversion relation for Li3. After subtracting the terms we 
have found in Eq. (6.4), we look at the image of the difference under A2,i or Ai^. As an 
example, we obtain 



A 



1,2 



U; ( - 

X 



T • / \ 1, ■! in , 2 

Li 3 (xH — In x m x 

V ' 6 2 



In 1 



In 2 



+ \ Mi 



x) (X> 



In x 



1, l2 * 7r 

— m x ® m x H 

2 2 



In x 



(6.5) 



0. 



We see that acting with A^2 does not provide any new information. This is not surprising, 
as the missing terms are of the form 7r 2 lnx, and Ai 5 2(7r 2 lnx) = 0. Indeed, acting with 
A2 1 and using the inversion relation for Li2, we obtain new non-trivial information, 



A 



2.1 



Li 3 



ITT 



Li 3 (x) H — In x In x 



Li 2 - 

x 



In I - 

X 



6 



Li2 (x) ® In x — - In x ® In x + i (it In x) ® In x 



1 vr 2 
Li2 (x) In x + inlnx -\ 

2 3 

1 



lnx 



(6.6) 



L12 (x) ® In x In x ® In x + i (it In x) ® In x 



1 



— 7T ® In x 
2 



Thus, 



3 

A 2 ,i 

Li.q 



7T 



lnx ) . 



1 „ 
Li 3 (x) + -ln d x 
6 



ITT 



In x 



7T- 



ln x + 0^3 + (3 m . 



(6.7) 



Specializing to x = 1 gives a = /3 = 0, which is indeed the correct inversion relation for 
Li 3 . Proceeding in exactly the same way, we can now derive the inversion relations for all 
the classical polylogarithms. 

6.2 Special values in x = 1/2 

As a second example we consider the special values of some harmonic polylogarithms when 
the argument is equal to 1/2. In many cases these values are expressible through £ values, 
In 2 and Li n (|), for n > 4. It is however impossible to obtain these relations using symbols 
alone, because 



H I 01, . . . ,a n ; 



-l) p 2®...®2 = 5 



(-iy 



iv. 



In™ 2 



(6.8) 



where ai £ {0, 1} and p is equal to the number of cij's equal to zero. As a consequence, a 
pure symbol approach only provides trivial and misleading information, because we always 
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obtain a symbol corresponding to powers of In 2. In the following we show that using the 
coproduct approach we can do better and entirely fix the values in x = 1/2, up to primitive 
elements of a given weight n (in the present case we only need to consider ( n ). 
We again proceed in a bootstrap and start from weight 2. We obtain 



Ai, 



Li 2 



In 1 



1)^(1 



■In 2(g) In 2 



Ax,! (In 2 2). (6.9) 



Thus, 



Li 2 



-In 2 2 + C7r 2 



(6.10) 



for some rational number c that cannot be fixed from the coproduct. Hence, at this stage 
we need to resort to numerics, 



Li 2 ( - ) + - In 2 2 = 0.82246703342411321824 ... = — 
'2/2 12 



(6.11) 



It appears at this stage that we have not gained anything over the pure symbol approach. 
In fact, in this case the coproduct only becomes more powerful than the pure symbol 
approach starting from weight three. Acting with Ai^i on L13 (2) yields 



Ai,i, 



Li 3 



In 2 <g> In 2 <g> In 2 = Ai 1 1 



In 3 2 



(6.12) 



Next, we can look at the (1,2) and (2,1) components of the coproduct, 



Ai, 



A 



2.1 



Li 3 


3) 


--ln 3 2 
6 


= 0, 


Li 3 


3) 


-iln 3 2 
6 


= -Li 2 ( - j ® In 2 - - In 2 2 <g> In 2 








= -— 7r 2 <g>ln2 = A 2 ,i 


In 2 

12 



Thus, we obtain 



Li 3 



Il n 3 2 _!Lln2 + aC3 + ^vr 3 . 
6 12 



(6.13) 



(6.14) 



j3 is obviously zero, because L13 (^) is real. Furthermore, from numerics we obtain 



U 3 [- 



Ii n 3 2 _!Li n 2 

6 12 



7 
1.0517997902646449997 . . . = - Cs ■ (6-15) 



We could now be tempted to apply the same procedure to Li4 Q) and express it in 
the form 

Li 4 ( I ) = a In 4 2 + c 2 vr 2 In 2 2 + c 3 Cs In 2 + c 4 n 4 . (6.16) 
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However, no such formula is currently known, and it is commonly believed that, starting 
from n = 4, Li n (2) defines a genuinely new transcendental number. If our 'extended 
symbol calculus' is consistent, it should lead us to the same conclusion, i.e., that an ansatz 
of the form (6.16) is excluded. To see why this is indeed the case, we start from Eq. (5.2) 
and we write 

U £ iZiL Li^* ( £ ) ® In* 2 . (6.17) 



A' 



Li 



a 



k\ '"12 
fc=i v 



We see that the second factor in the reduced coproduct only involves powers of In 2. An 
ansatz made out of combinations of powers of In 2 and £ values however inevitably leads 
to terms in the coproduct that have a C, value in the second factor, e.g., 

A(C m ln fe 2) = A(C m ) A(ln fc 2) = ( m ln fc 2 + ln fe 2 <g> ( m + . . . . (6.18) 

The only way to make the terms having a ( value in the second factor vanish is to assume 
that m is even, because Eq. (5.15) implies 

A(^ m ln fc 2) =7r m ®ln fc 2 + ... . (6.19) 

Thus, we conclude that any ansatz made out of products of powers of In 2 and £ values 
can only involve even £ values. This is indeed what happens for Li3 (2) where the only 
possibility for a product of the form Cm l n 2 is 7r 2 ln2. Starting from weight four this is 
excluded, because Eq. (6.17) involves a term 

Li 3 (^) In™- 3 2 = I Ca ® hi™" 3 2 + . . . . (6.20) 

We therefore arrive at the conclusion that starting from weight four, Li n (^) can no longer 
be expressed through ( values and powers of In 2 alone, in agreement with the common 
belief. We stress the role played in this argument by the special treatment of £2, Eq. (5.12). 
So far we have only considered examples of classical polylogarithms. Let us therefore 
conclude this section by discussing a less trivial example of weight five, where the full 
superiority of the coproduct approach over the pure symbol approach is revealed. Consider 
to this effect the harmonic polylogarithm H (0, 1,0,0, 1; 2)- We can make an ansatz for 
this number in the form 

T =ciLi 5 ( - J +c 2 Li 4 ( - J ln2 + c 3 In 5 2 



2/ W (6.21) 

+ c 4 it 2 In 3 2 + C5C3 In 2 2 + c 6 7T 4 In 2 + c 7 vr 2 (3 + c 8 C5 • 

Our goal is to find rational numbers q such that H (0, 1, 0, 0, 1; i) = T. A pure symbol 
approach is obviously totally inadequate for this: not only would it be unable to constrain 
the coefficients multiplying the £ values, but it could also not distinguish between the first 
three terms in the ansatz, thus only providing a single relation between the coefficients c±, 
C2 and C3. As we will see in the following, the coproduct approach allows us to fix all the 
coefficients, except for the coefficient of £5 (which is a primitive element). 
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We start by computing the maximal iteration of the coproduct, 

Ai,i, 1,1,1 (T) = (-ci + 5c 2 - 120c 3 -I)ln2(8)ln2(8)ln2(8)ln2(8)ln2, (6.22) 

where we introduced the shorthand 

T = .fffo,l,0,0,l;^-T. (6.23) 

Equating the right-hand side to zero, we obtain a relation between c±, c 2 and C3, 

Cl = 5c 2 - 120c 3 - 1 . (6.24) 

We note that this relation is precisely the information we could extract from a pure symbol 
approach. We can however now go on and compute the components of the coproduct 
that involve precisely one factor of weight w > 1. Inserting the solution (6.24) into the 
expression for T, we obtain 



and 



and so we can solve for c 2 , C3 and C5, 



For the (3,2) and (2,3) components of the coproduct we obtain 



.25) 



Ai.i.1,2 (T) = Ai,i, 2 ,i (T) = Ai, 2 ,i,i (T) = 0, 

A 2 ,i,i,i (T) = (— - 10c 3 - 6c 4 ) vr 2 <g> In 2 ® In 2 ® In 2 , 

and so 

c 4 = — - -c 3 • (6.26) 

36 3 

Next we investigate the components of the (iterated) coproduct m.'H§$>'H® / H, 

Ai, 2 ,2 (T) = A 2 ,i, 2 (T) = A 2i2 ,i (T) = , (6.27) 



A; L .j., (V) - ( J-2c 5 jln2^1n2®C3, 

Ai.^lT) - (^-2c 5 -^ln2®C3®ln2, (6.28) 

A:u.j (■'/') - (-^ + 105c 3 -2 C5 V3®ln2®ln2, 



11 7 

C2 = 3 and C3 = and cs = - . (6.29) 

120 8 



A 3 , 2 (T) = and A 2;3 (T) = - I c 7 + — J tt 2 (8> Cs , (6.30) 

and so C7 = — jg. Similarly the (4,1) and (1,4) components yield 

Ai, 4 (T) = and A 4 ,i (T) = f ^g " c 6 J vr 4 In 2 , (6.31) 
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and so c% = ^gg- Finally we arrive at 

H (o, 1,0,0, l;i) -SU.(i) + Sh.2U.(i)+^h»2-^W2 

7 1 7 l j 
+ - C 3 In 2 2 + vr 4 In 2 vr 2 C3 + c 8 (5 ■ 

8 s 288 48 s s 

As expected, the coproduct allowed us to fix all the coefficients except for c 8 . Using 
numerics, we arrive at 

H (0, 1, 0, 0, 1; ^ - T = -eg C 5 - 1.3123616901033275 ...=-<* Cs - j| Cs , (6-33) 

and thus c 8 = — |^. 

7. Amplitudes for H + 3 gluons 

In this section we apply the coproduct to a physical problem, namely the two-loop helicity 
amplitudes for a Higgs boson plus three gluons in the large top mass limit. In this limit 
the coupling of a Higgs boson to gluons is described by an effective operator of dimension 
five, 

C eff = -^HG%G^. (7.1) 

The two-loop corrections to the helicity amplitudes for a Higgs boson plus three gluons were 
computed in Refs. [60, 61], where it was expressed as a complicated combination of two- 
dimensional harmonic polylogarithms. In Ref. [53] it was shown that, after subtracting the 
square of the one- loop amplitude, the symbol of the leading color maximally transcendental 
part of the two-loop helicity amplitudes is equal to the symbol of the two-loop form factor 
of three gluons in planar J\f = 4 Super Yang-Mills. The latter can be expressed in a very 
compact form involving only classical polylogarithms up to weight four [53]. This suggests 
that the two-loop corrections to the helicity amplitudes for a Higgs boson plus three gluons 
can be written in a much simpler form without any multiple polylogarithms. However, as 
the symbol does not fix terms proportional to £ values, the symbol alone is insufficient to 
determine such a simplified form in an easy way. In the following we apply our coproduct 
approach to rewrite the results of Refs. [60, 61] in a compact form, obtaining in this way 
compact analytical expressions for all helicity amplitudes for a Higgs boson plus three 
gluons, for both the decay {H — > ggg) and the scattering (gg — >■ Hg) regions. 

7.1 The decay region 

We start by investigating the decay region, i.e., the two-loop corrections to the helicity 
amplitudes for H — > ggg. The kinematics is described by three dimensionless ratios, 

•512 S 2 3 S31 . . 

xi = —5- , x 2 = -^r , %3 = -^r , (7.2) 

m H m H m H 

where run denotes the mass of the Higgs boson and Sij = 2piPj , with pi the momenta of the 
external gluons. These kinematic variables are not independent, but they are constraint 
by 

< Xi < 1 and X\ + X2 + X3 = 1 . (7.3) 
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As a consequence, the amplitude is effectively a function of only two of the three dimension- 
less ratios. Correspondingly, the result of Ref. [61] is expressed in terms of two-dimensional 
harmonic polylogarithms in x 2 and x 3 . There are two independent helicity configurations 
for the decay, 

H^g + g + g + and H -> g + g + g~ . (7.4) 

In the following we will analyze each configuration separately. 

Let us start by analyzing the helicity amplitude where all the final state gluons have 
a positive helicity. Bose symmetry then implies that the amplitude must be symmetric 
under a permutation of the external gluons, or, equivalently, it must be totally symmetric 
in the kinematic variables Xi, i = 1,2, 3. The (finite part of the) one-loop correction to the 
decay can be written as 

a™ = N C A^ + (^ - H Nf \ B& + \{N C - N f )( Xl x 2 + x 2 x 3 + x 3 x x ) , (7.5) 



with 



A^> = — -(lnxi lnx 2 + hix 2 lnx 3 + lnx 3 lnxi) - y^Li 2 (l - xi) , 

4 2 7^1 ( 7 - 6 ) 

B^> = ln(xi x 2 x 3 ) - Zisi . 

Following Ref. [61], we decompose the (finite part of the) two-loop correction into contri- 
butions with different color structures, and we furthermore subtract the square of the finite 
part of the one-loop amplitude, Eq. (7.5), 



«( 2 ) = - r 



a« 



+ Nl if + ^ Jjf + Nc Nf E^ + N] Fl 2) . (7.7) 



The coefficients of the different color structures were computed in Ref. [61] where they 
were expressed as a combination of two-dimensional harmonic polylogarithms in x 2 and 
X3. In order to simplify these expressions, we start by computing the symbol of Eq. (7.7). 
It turns out that all the entries in the symbol are drawn from the set 

{xi, x 2 ,x 3 , 1 - xi, 1 - x 2 , 1 - x 3 } . (7.8) 

The weight four part of the symbol satisfies 

" =0, (7.9) 



S ( a (2) 



weight 4 

where 5 (a <8> b <S> c <8> d) = (a A b) A (c A d), and the wedge denotes the antisymmetric tensor 
product, a /\b = a®b — b<S> a. It follows then from a conjecture in Ref. [66] that cr 2 ' can 
be expressed in terms of classical polylogarithms only. Similar conclusions were already 
drawn in Ref. [53]. 

Next, we have to determine the arguments of the polylogarithms that can lead to 
a symbol with entries drawn from the set (7.8) under the constraint (7.3). Using the 
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XI 


1 — Xl 




1 - 1/xi 




X2 


1 - X 2 




1 - l/x 2 




x 3 


1 -x 3 




1 - l/x 3 




-xi/x 2 


^2/(1 -x 3 ) 




Xi/(1 -x 3 ) 




~X2/X 3 


x 3 /(l -xi) 




£2/(1 -xi) 




-X 3 /X! 


»i/(l -x 2 ) 




£3/(1 -Z2) 




~X 1 X 2 /x 3 


x 3 /[(l-xi)(l- 


*2>] 


XlX 2 /[(l -xi)(l - 


■X2)} 


-X 2 X 3 /xi 


Xi/[(l-x 2 )(l- 


*s)] 


Z 2 X 3 /[(1 - X 2 )(l - 


■xs)} 


-x 3 xi/x 2 


x 2 /[(l-x 3 )(l- 


*l)] 


x 3 xi/[(l -x 3 )(l - 


■xi)} 



Table 1: Arguments of classical poly logarithms that can give rise to a symbol with entries drawn 
from the set in Eq. (7.8) under the constraint (7.3). Each line shows half an orbit of the S3 action, 
the second half being obtained by inversion. All these functions are less than unity in the region 
defined by Eq. (7.3). 



prescription given in Ref. [57], we find 54 rational functions grouping into 9 orbits of the 
symmetric group 5 3 whose action on rational functions /(xi,x 2 ,x 3 ) is generated by 5 



/->!-/ and /-►!// 



(7.10) 



The rational functions are summarized in Table 1. It is important to note that not all 54 
solutions are independent, and in particular we can express half of them in terms of the 
others by using the inversion relation for the classical polylogarithms, 

1 



Li T 



ir +i Li„(/) + .... 



(7.11) 



It is then easy to see that it is always possible to choose 27 solutions such that all polylog- 
arithms are real in the region denned by Eq. (7.3). 

Next step we write down a combination of (classical) polylogarithms in the arguments 

shown in Table 1. Equating the symbol of cr 2 ) and our ansatz provides a linear system for 

(2) 

the coefficients. In the following we only discuss the weight four part of A a . All other 

contributions are similar. In agreement with Ref. [53], we find 



S A 



r(2) 

a, weight 4 



sin 



,(2) 



(7.12) 



,(2) 



where 7^.3 is the M = 4 form factor remainder function of Ref. [53] , 

3 



n 



(2) 



1 

12 



A, 



XlX 2 
^3 



.4 = 1 v ' 



+ A 4 
2 



£1^3 
X2 



+ A 4 



X2X 3 \ 
Xi / 






2 1 

-lnxi lnx 2 lnx 3 ln(xix 2 x 3 ) H In (xix 2 x 3 ) 

3 16 



(7.13) 



1 ^^TT^ 

H — (lnxi lnx 2 + lnxi lnx 3 + lnx 2 lnx 3 ) (in xi + In x 2 + In x 3 ) — , 



3 We stress that this S3 symmetry is not identical to the S3 describing the Bose symmetry. 
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where A n (z) denotes Rummer's function 

In"" 1 \t 



K(z)=l Q dt- ]+t 



" 1 ' i \n~ k 



("~1)'E m ln fe klLin-fe(^)- 



fc=0 



A:! 



(7.14) 



This result was already obtained in Ref. [53]. However, Eq. (7.12) only holds at the level 
of the symbol, and it would thus be premature to conclude that the weight four part of 

(2) (2) 

A a is equal (at the level of the function) to TZ 3 . Indeed, acting with A2,i,i, we obtain 



A 2 ,i,i 



A 



(2) 

a, weight 4 



Tt 



(2) 



7T <g> Ai i 



A W 



A 2 ,l,l 



7T~ 



■*■ - L /-v 



Continuing this way, we can easily determine the coefficient of £3, 



A 



3,1 



-T(2) 

a, weight 4 



1 



*?> + y 4 X) ] = 4 Ca ® i#> = A 3 , [ -| Ca < 



(7.15) 



(7.16) 



Finally, we determine the coefficient of 7r 4 by evaluating the function at a single point in 
phase space, 



r(2) 



71" 



1 



-^•a, weight 4 3 fi a A ^ 3 a 



-0.03382260105347.. 



7T 



2880 



(7.17) 



Repeating the same steps for all other contributions to Eq. (7.7), we arrive at the following 
expressions for the different color structures contributing to the two-loop amplitude aS 2 > , 
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a ~^ 3 6 a 4 C3 a 2880 
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+ - ln(xi x 2 x 3 ) A^ + - ^2 [Li 2 (1 - Sj) In Xi] + - In xi In x 2 In x 3 + - In 3 (xi 



S2S3J 
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x? iX 2 
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P 2 (Xi,Xi-i,X i+ i) 



121 1 2 
mxi_i mxi + i H — In x. 



144x 2 x 2 x 2 
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216~' 



(7.18) 
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£>a = -C3 + T ~ fi ( XlX2 + X3X2 + XlXZ > + ~ + — '"' " '"' " 
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Pg(Xj,Xj_i,Xj + i) 
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(7.20) 
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18 v ; 144 36 4^ 54 v ; 18 



i=l 



m , . 5 . . 
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72 



(X1X2 + X3X2 + X1X3 



2 Xix 2 x 3 ^— v lnx,- 



18 



i=l 



(7.21) 



where Pj(x, y,z) = Pi(x,z,y) are homogeneous polynomials in three variables, 



Pi(x, y, z) = 30x 2 (y 4 + z 4 ) - 134y 2 z 2 (x 2 + y 2 + z 2 ) - 199xy 2 z 2 (y + z) 

+ 75xyz (xy 2 + xz 2 + y 3 + z 3 ) - 268y 3 z 3 , 
P 2 (x, y, z) = — 134x 3 (y + z) — 67x 2 (x 2 + y 2 + z 2 ) — 65x 2 yz + 75xyz(y - 
P%{x, y, z) = —20 (x y + x z + y z ) — 276x y 2 (xy + xz + yz) 

/ 9 9 9 9 9 9 99 99 9 9\ 



= -20 (x 4 y 

••'' - 2 ' x 2 y 3 + x 2 z 3 + y 3 z 2 + y 2 z 



- 50xyz (x 3 y 2 + x 3 z 2 + x 2 y 3 

- 115x 2 y 2 z 2 (x 2 + y 2 + .- 2i 



,2^3\ 



P 4 ( x , y, z) = 60x 2 (y 2 + z 2 ) + 9x 2 yz - 78xyz(y + z) - 256y 2 z 2 - 117yz (y 2 + z 2 ) , 
P 5 (x, y, z) = x 4 y 4 + x 4 z 4 + y 4 z 4 - 2xyz (x 3 y 2 + x 3 z 2 + x 2 y 3 + x 2 z 3 + y 3 z 2 + y 2 z 3 ) , 



(7.22) 
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P 6 (x, y, z) = -x [x (y A + z A ) - 2xyz (y 2 + z 2 ) - 2yz (y 3 + z 3 )] , 

Pl{x, y, z) = yz(2xy + 2xz - yz) , 

P 8 (x, y, z) = -x [2x (y 2 + z 2 ) - hxyz - 4yz(y + z)] , 

P 9 {x, y, z) = -33x 2 {y A + z A ) + 20yV (x 2 + y 2 + z 2 ) - 2xy 2 z 2 (y + z) 



A2xyz (xy 2 + xz 2 + y 3 + z 3 ) + 40y 3 z 3 , 



Pio(x, y, z) = 20x 3 (y + z) + 10x 2 (x 2 + y 2 + z 2 ) - 22x 2 yz - 42xyz(y + z) - 33y 2 z 2 , 
P\\ (x,y,z) =44 (x A y A + x A z A + y A z A ) + 410x 2 y 2 z 2 (xy + xz + yz) 



+ 56xyz (x 3 y 2 + x 3 z 2 + x 2 y 3 + x 2 z 3 + y 3 z 2 + y 2 z 3 ) 

+ 177x 2 y 2 z 2 (x 2 + y 2 + z 2 ) , 
Pu(x, y, z) = -198x 2 (y 2 + z 2 ) + 35x 2 yz + U2xyz(y + z) + 490y V + 206yz (y 2 + z 2 ) , 
Pi 3 (x, y, z) = -1781 (x 4 + y A + z A ) - 8224 (x 3 y + x 3 z + xy 3 + xz 3 + y 3 z + yz 3 ) 

- 12874 (x V + x 2 z 2 + y 2 z 2 ) - 26848xyz . 

(7.23) 

Note that P3, P5, Pn and P13 are totally symmetric. The expressions for the two-loop 
corrections to cr 2 ' in Eqs. (7.18 - 7.21) only involve classical polylogarithms with a rather 
simple dependence on the kinematic invariants. In particular, the functional dependence 
is such that all the polylogarithms are real in the region (7.3). Furthermore, the S3 Bose 
symmetry of the amplitude if completely manifest in our expressions. 

We can of course apply exactly the same procedure to the second helicity configuration, 
H —?■ g + g + g~ . Bose symmetry implies that this amplitude must be symmetric under an 
exchange of the two positive-helicity gluons, or, equivalently, under an exchange of X2 and 
X3. The one-loop corrections are given by 



/3« = N c Af + (^--N f 
H c P \ 6 12 ; 



with 



,(1) 

A f3 



*?' +!<".- 


X\ 


(7.24) 


3<" = B m . 




(7.25) 



For the two-loop corrections we write 



/5 (2) = \ b (1) l 2 + N 2 C Af + ^ Df + N c N f pf + N 2 pf , (7.26) 



2L r J ' " u "" p ' iV L 
Repeating the same steps as for the first helicity configuration, we arrive at 

/3 3 g a 4 S3 a 2gg0 g a t gQ 



+ yLi 3 (l-x) + yLi 3 I 1--J - ^gg*" 2 Inxi -— ln(x 2 x 3 )Li 2 (l-xi) 

- — In 3 xi + — - In 2 xi ln(x 2 x 3 ) + — In 2 x 2 hi(xix 3 ) + — In 2 x 3 ln(x 2 xi) 
6b 24 24 24 



(7.27) 
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Li 2 (1 - xi) + 



Q25(x 1 ,X 2l X 3 ) 2 



79„,4„,2„,2 



-IT 



(l-x 2 ) x 2 Q 24 (xi,x 2 ,x 3 ) (I-X3) x 3 Q 24 (xi,x 3 ,x 2 ) T . ,., , 

j-^ Ll 2 (1 - X 2 ) 1n -4-2 Ll 2 I 1 - X 3j 



Ylx\x\ 



Ylx\x\ 



1 

12 



lnxi 



x 2 (x 2 - 2x 3 ) x 3 (x 3 - 2x 2 ) 
5 ln x 2 H 5 ln x 3 



+ 



Q 22 (xi,x 2 ,x 3 ) 
48x 2 (l-x 2 )(l-x 3 ) 

x 2 x 3 (6xf + 4x 2 xi + 4x 3 xi + 3x 2 x 3 ) (2x 2 - x 3 )(2x 3 - x 2 ) 
2 m x 2 In x 3 -\ — in x 4 



+ 



12xf 
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, x 3 Q 23 {xi,x 3 ,x 2 ) 
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F (2) - -™ ln^xox^ + — Vln 2 x- - — vr 2 5 (^2 + 4) + 13a ^ 2 - 10(x 2 + x 3 ) 

*> - i8 ln(xix 2 x 3 ) + 36 2^1n x, ^tt lnxi 



+ 



5(x 2 + x|) + llx 3 x 2 - 10(x 2 + x 3 ) + 5 



18x 2 



X 2 X 3 / 2 \ 5 

"r + — -j ( 20x i - 3x 2 x 3 ) - — ln(xix 2 x 3 ) 

(7.31) 
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where the Qi(x,y,z) are homogeneous polynomials, 
Qi(x,y,z) ■■ 



39x 5 z + 20x V + 48x 4 yz - 78x 4 z 2 + 30x 3 y 2 z - 20x 3 yz 2 
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Q22 


[x 


y, 


z) 
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2 A 
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24xy 2 z 3 - 24xyz 4 + 90 y z z 4 , 



30y 4 + 75y 3 z - 134y 2 z 2 + 75yz 3 + 30z 4 . 
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50x 4 y z 3 



20xV - 4x 3 y 3 z 2 



4x 2 y 2 z 4 + 16xy 4 z 3 + 16xyV 
206x 5 z - 198xV - 234x 4 yz + 412x 4 z 2 - 297 x 3 y 2 z + 72x 3 yz 2 



60yV 



+ 206xV 



I26x 2 y 2 z 2 + 342x 2 yz 3 



855 xy 2 z 3 + 36xyz 4 - 594y V 
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-6x 3 



3x 2 y 



3x z 
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2y 3 , 
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6x 2 y 2 z 2 + 48x 2 yz 3 



+ 6xV + 12xy 2 z 3 + 12xy z 4 



99y 2 z 4 , 



10x 4 + 48x 2 yz + 12xy 2 z + \2xyz 2 



99y 2 z 2 



= 44xV + 56x 4 y 3 z + 177x 4 y 2 z 2 + 56x 4 y z 3 + 44xV + 8x 3 y 3 z 2 + 8x 3 y 2 z 3 
+ 8x 2 y 4 z 2 - 64x 2 y 3 z 3 + 8x 2 y 2 z 4 - I6xy 4 z 3 - 16xy 3 z 4 + 132y V , 



26xy z — 26xyz 



67x 4 + 65x 3 y + 65x 3 z - 2x 2 y 2 + 27x 2 yz - 2 x 2 z 2 

\2y 2 z 2 , 

2x 4 y - 4x 4 z + 3x 3 yz + \2x 3 z 2 + \8x 2 y z 2 + 24x 2 z 3 + 17xyz 3 + 8xz 4 



+ 6yz , 



x 2 y 



2x 2 z — 2xyz + 4xz 2 + 3yz 2 , 



x 4 y 4 - 2x 4 y 3 z - 2x 4 yz 3 + x 4 z 4 + 6x V z 3 + 4xy 4 z 3 + 4xy 3 z 4 + 3y 4 z 4 . 

(7.32) 

As expected, the results involve only classical polylogarithms with arguments such that all 
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functions are real in the region (7.3). Furthermore, the Z2 Bose symmetry is completely 
manifest. However, due to the reduced symmetry with respect to the first helicity configu- 
ration, the expressions are not as compact as in the previous case. It is worth noting that 
the weight four contribution is identical for both helicity configurations. 

7.2 Analytic continuation to the scattering region 

The expressions presented in the previous section are only valid in the decay region (7.3). In 
the rest of this section we show how to perform the analytic continuation to the scattering 
region. We have to distinguish the following cases, 

g + g + -> Hg + } 

TT > xi > and x 2 ,x 3 < 0, 
g + g + ^Hg J (7.33) 

g + g~ — > Hg + x 2 > and X\,x 3 < 0. 

In all cases the kinematic invariants are subject to the constraint 

xi + x 2 + x 3 = 1 , (7.34) 

which simply expresses s + t + u = m 2 H . In the following we only discuss the analytic 
continuation in the case where all gluons have a positive helicity, all other cases being 
similar. 

In the decay region all invariants are positive and have a small positive imaginary 
part. The analytic continuation to the scattering region is then performed according to 
the prescription 

s 2 3 -> |«2s| e in and s 13 ->• |s 13 | e ilT , (7.35) 

while all other invariants remain unchanged. This implies the following prescription for 
the dimensionless ratios, 

x\ — > x\ and x 2 — > x 2 e l7T and X3 — > x 3 e ln , (7.36) 

where we defined Xj = |xj| = — X{. 

Using these prescriptions, the Kummer functions are analytically continued according 
to 

A„(-ze*) -)• (n - 1)! ^ [ ~ ^ [In \z\ + iJvr] fc Li TO _ fc (-z) . (7.37) 

fc=0 
In addition, we need to analytically continue classical polylogarithms of the form 

Li n (l-z e l& A , O0 and <5 = ±1. (7.38) 

While the corresponding formulas could be obtained by the help of, e.g., the Mathematica 
package HPL [32], we show how the analytic continuation formulas can be derived from 
the coproduct. Similar to the case of the inversion relations discussed in Section 6, we 
proceed recursively in the weight. At weight one, we immediately obtain 

Lii (l-z e i57T ) = - In (z e iS A = - In z - iba . (7.39) 
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At weight 2, we act with the coproduct, and drop all the in terms in all the factors of the 
coproduct except the first one, 



A 



1.1 



Li 2 l-ze l 



Lii 1 - z e tdw <g> In 1 - z e 



ASn 



In z <g> ln(l + z) - in <g> ln(l + z) 



(7.40) 



A i,i 



Li 2 (— z) — lnzln(l + z) — Mrln(l + z) 



Thus, we obtain 



Lio ( 1 — z e 



'■()7T 



-Li2(— z) — lnzln(l + z) — mln(l + z) + en , 



(7.41) 



for some rational number c. Specializing to z = (where we are insensitive to the phase) 
we immediately obtain c = g. 

At weight 3, we first act with Am, 
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Li, (1-ze* 



ASt. 



iSn 
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,i57 
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1,1,1 



Lia 



1 + z 



-ln 3 (l + z)--(57rln 2 (l + z) 
6 v 7 2 v 7 



(7.42) 



In order to determine the terms proportional to n 2 , we compute 
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Li 3 (l-ze^)- 

= -n 2 ® ln(l + z) 
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= A 2 ,i 


M x iv 


1 + zJ 6" 

— ln(l + z) 



^ln 3 (l + z)--57rln 2 (l + z; 



(7.43) 



Thus, 



Li 3 1 - z e 



Li 3 



1 + z 



1 i n 2 

-ln 3 (l + z)--5vrln 2 (l + z) + — ln(l + z) + cC 3 . (7.44) 
o 2 3 



Specializing to z = 0, we obtain c = 0. 

Using this technique we can recursively derive all the analytic continuation formulas 
for functions of the type (7.38). In particular, at weight 4 we obtain 



Li 4 I 1 - z e l 



-Li/i 



1 



1 



i5n 



n- 



n 



ln 4 (l + z)-^^ln 3 (l + z) + — ln 2 (l + z) + — . (7.45) 
1 + z/ 24 v ' 6 v ' 6 V ' 45 v ; 



These formulas are enough to perform the analytic continuation from the decay region 
to the scattering region. We checked numerically that our results agree (after analytic 
continuation) with the results in the scattering region presented in Ref. [61] for all the 
helicity configurations. 
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8. Conclusion 

While recent advances seem to indicate that, at least in the context of the M = 4 Super 
Yang-Mills theory, scattering amplitudes are simpler than expected, it is a known fact 
that the analytical evaluation of multi-loop Feynman integrals can lead to very lengthy 
and complicated expressions involving new classes of transcendental functions only poorly 
studied in the literature. A systematic approach to study these new functions and their 
functional equations is therefore highly desirable, not only from the formal standpoint, but 
also in perspective of phenomenological applications. A first step in this direction has been 
made in Ref. [40] with the introduction of the so-called symbol map that allows to map the 
combinatorics of transcendental functions defined by iterated integrals to the combinatorics 
in a certain tensor algebra. 

In this paper we have proposed a novel approach to deal with complicated expressions 
that can arise from a special class of Feynman integral computations, namely those that 
can be evaluated in terms of multiple polylogarithms. The cornerstone of this approach is 
the coproduct on multiple polylogarithms introduced by Goncharov in Ref. [58], augmented 
by some ideas from a recent paper by Brown [59]. The main feature is that, unlike the 
symbol, the coproduct allows to incorporate £ values into the calculus, thus retaining more 
information about the function. We have demonstrated the virtue of this novel approach 
by rewriting the two-loop helicity amplitudes for a Higgs boson plus three gluons originally 
computed in Ref. [60, 61] in a compact analytical form, revealing, at least to our knowledge, 
for the first time an unexpected simplicity for a two-loop multi-scale amplitude in QCD. 
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A. The multiple polylogarithm coproduct 

A.l The coproduct in the generic case 

Before discussing how to define the coproduct of a multiple polylogarithm with non-generic 
arguments, let us first review a diagrammatic interpretation of the coproduct introduced 
in Refs. [58, 63]. 

Consider a multiple polylogarithm of the form I{clq; a\, . . . , a n \ a n+ \). We associate to 
it a diagram where we arrange the arguments of the multiple polylogarithm on a semi-circle. 
For example, for n = 3 we have 




(A.l) 
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The different terms in the coproduct in Eq. (5.1) then correspond to connecting points 
via a polygon (including the empty polygon) in all possible ways. The points lying on the 
polygon provide the arguments for polylogarithm in the first factor in a given term in the 
coproduct, while the remaining points determine the entry of the second factor. Here we 
illustrate this construction only on the example of I(ao; 01,02,03, 04; 05), and we refer to 
Refs. [58, 63] for further details. In this case the polygons together with the terms in the 
coproduct they correspond to are 




a a 5 

1 <g> I(a ; ai,a 2 , a 3 , a 4! 05) 




a a 5 

I(a ; ai; 05) <g> I(oi; o 2 , 03, a 4 ; a 5 ) 



a 2 03 



Cl2 0,3 

a\S y/ \\ a4 



ao «5 a o 0,5 

I(a ;a 2 ;a 5 ) <g> [I (a®; ai;a 2 )I(a 2 ; 03,04; 05)] /(a ;a 3 ;a5) [/(oo; 01,02; 03)1(03; 04; 05)] 




a a 5 

7(a ; 04; as) <8> /(a ; ai,a 2 , 03; 04) 




ao 05 

I(a ; ai,a 2 ; a 5 ) <g> I(a 2 ; a 3 , a 4 ; a 5 ) 



a 2 a 3 

ai f^^^ \\ a 4 

r ^ 

a a 5 

I(a ; ai,a 3 ; a 5 ) <g> [J(ai; a 2 ; 03)^(^3; a 4 ; a 5 )] 




ao 05 

I(a ; ai,a 4 ; a 5 ) (g) J(ai; a 2 , a 3 ; a 4 ) 
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J(a ; 02, 03; «5) <S [/(a ; an a 2 )I(a 3 ; a 4 ; a 5 )] I(a ; a 2 , a 4 ; a 5 ) [I(a ; an a 2 )I(a 2 ; a 3 ; a 4 )] 



a 2 03 

ai <r ^/ >, a 4 

/ \ 

a a 5 

I(ao;a3,a4;a 5 ) <g> I(a ; ai, a 2 ; a 3 ) 




a as 

7(a ; 01,02, 03; as) (S> /(a 3 ; a 4 ; 05) 




a a 5 

7(a ; 01, a 2 , a 4 ; a 5 ) <g> I(a 2 ; a 3 ; a 4 ) 



a 2 03 

r a 

7(a ; ai,a 3 , a 4 ; 05) (S> /(an a 2 ; 03) 




00 05 

7(a ; 02, 03, «45 a 5 ) (2) I(a ; an a 2 ) 




a 05 

I(a ;ai,a2,a 3 ,a 4 ;a 5 ) <g> 1 



A. 2 The coproduct in the non-generic case 

As already mentioned, the formula for the coproduct, Eq. (5.1), is only valid in the 
generic case where all the arguments are mutually different. Indeed, in the non-generic 
case divergences can arise in individual terms in the coproduct. A multiple polylogarithm 
/(ao; ai, . . . , a n ; a n+ i) is in general divergent if either ai = ao or a n = a n+ ±. In this case the 
poles in the integrand coincide with the endpoints of the integration path. As an example 
of how these divergences can arise inside the coproduct, consider the multiple polyloga- 
rithm /(ao; ai, a 2 , a 2 ; a 3 ), which is convergent whenever ao 7^ ai and a 2 / a 3 . Eq. (5.1), 
however, contains a term 

a 2 




a a 3 

J(a ;ai,a 2 ;a 3 ) <g> /(an a 2 ; a 2 ) . 
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and the second factor, 7(ai;a 2 ;a 2 ), is divergent. 

In Refs. [58, 63] the coproduct in the non-generic case is defined by replacing in the 
right-hand side of Eq. (5.1) every multiple polylogarithm by its regularized value. In the 
following we only give a practical rule of how to obtain a regularized value (there is more 
than one way to perform the regularization) , and we refer to Refs. [58, 63] for more details. 

As the divergences of a multiple polylogarithm are end-point divergences, we can easily 
regularize them by slightly moving the end points of the integration path. In practice we 
usually deal with integration paths that are straight lines, and in this case the regularization 
can simply be achieved by the replacement 



I(ao;ai, • • • ,a n ;a n +i) 



I(a (l + e); oi, . . . , a n ; a n+1 (l - e)) , if a / . 
I(e; ai , . . . , o n ; a„+i (1 - e)) , if a = , 



(A-2) 



This regularization procedure has two obvious features: 

1. If I (ao; a±, . . . , a n ; a n +i) is convergent, then 

lim I(oq(1 + e);ai,... ,a n ;o n+ i(l - e)) = J(a ;ai, . . . ,a n ;a n+ i) . 

2. The regularization (A. 2) preserves shuffle identities. 



(A-3) 



As a consequence, we can use shuffle identities to extract all the divergences of the multiple 
polylogarithms, e.g., 

I(a (l +e);a 1 ,a 2 ,a 2 ;a 2 (l - e)) = I(a (l + e);a 2 ,a 2 ;a 2 (l - e))I{a (l +e);oi;o 2 (l - e)) 
-I(oo(l + e);o2,oi,02;o 2 (l - e)) -I(a (l + e);a 2 ,a 2 ,ai;a 2 (l - e)) 
= I(oo(l + e);a 2 ,a 2 ;a 2 (l - e))I(a (l + e);ai;a 2 (l - e)) 
- I(a (l + e); a 2 , oi; o 2 (l - e))I(a (l + e); a 2 ; a 2 (l - e)) 
+ J(a (l + e); a 2 ,a 2 ; ai(l - e)) 



-ln 2 e/(a ;ai;a 2 ) - lne/(a ; a 2 , ai; a 2 ) + 7(a ; a 2 ,a 2 , oi; a 2 ) 



(A.4) 



In this way all the divergences are regularized, and the shuffle identities allow us to write the 
regularized integrals as a polynomials in lne. The regularized value /(ao; a±, . . . , a n ; a n +i) 
of the multiple polylogarithm /(ao; a±, . . . ,a n ; a n+ \) is defined as the constant term of this 
polynomial. As an example, from Eq. (A.4) we obtain 



J(ao;ai,a 2 ,a 2 ;a 2 ) = /(ao; a 2 ,a 2 , ai; a 2 ) . 



(A.5) 



The coproduct in the non-generic case is then obtained by replacing all multiple polyloga- 
rithms in the right-hand side of Eq. (5.1) by their regularized values, 



A(/(oo;oi,... ,a n ;a n+ i)) 

= 2_^ I{ a o) a iii ■ ■ ■ ^i k ;a n+ i) 

0=il<t2<---<ik<ik+i=n 



p=0 



+1) 



, a 



ip+i- 



-i;a 



lp+1, 



(A.6) 
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B. Proofs of the derivative and monodromy identities 

B.l Proof of the derivative identity 

In this section we sketch the proof of the derivative identity (5.23) in a particular case, 
namely 



id 



a 



da n+ i 



A(/(oo; 01, . . . , o n ; o n +i)) = A 







da 



n+l 



-I(a ;ai,...,a n ;a n +i) , (B.l 



where we assume all arguments of the multiple polylogarithm generic. 

Let us compute the action of the differential operator in the left-hand side of Eq. (B.l). 
It is obvious that we do only get a non-zero contribution from those terms in the coproduct 
where the second factor depends on a n+ \. If all the arguments are generic, it is easy to see 
that this is the case precisely for those terms in Eq. (5.1) where the polygon inscribed into 
the semi-circle (see Appendix A) does not contain a n . As an example, for n = 3, the term 




dQ <2 4 

I(a ;ai,a 3 ;a 4 ) <g) I(ai;a 2 ;a 3 ) . 



vanishes under the action of the differential operator ( id <g> 

a 2 




ao 04 

I(a ; ai; a 4 ) <g) I{a 1 ;a 2 , a 3 ; a 4 ) , 



gives a non-zero contribution, namely 



a 4 -a 3 




a 
, 7(a ;ai;o4 



/(oi;a 2 ;a 4 ) 



whereas the term 



(B.2) 



It is clear that in this way we produce precisely all the polygons inscribed into the semi-circle 

with the point a n (a 3 in the example above) removed, multiplied by = — . These terms 

are in one-to-one correspondence with the terms in the coproduct of /(ao; a±, . . . , a n -\; a n +i)- 
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Thus, we obtain 

d \ 1 

id ® a A(/(oo;oi,... ,a n ;a n+ i)) = A (I(a ; a±, . . . ,a n _i;a n+ i)) 

OO-n+l J «n+l — 0-n 

= A I I(a ; ai, . . . , a n _i; a„+i) ) (B.3) 

\ a n+l ~ a n J 

-I(ao;ai, . . . , a n ;a n+ i) , 



which finishes the proof. 

B.2 Proof of the monodromy identity 

In this section we sketch the proof of the monodromy identity (5.23) in a particular case, 
namely 

{Ma n+l = ai ® id) A(I(a ; ai,...,a n ; a n+1 )) = A (A4 an+1=ai /(a ; a±, . . . ,a n ; a n+1 )) , (B.4) 

where all arguments of the multiple polylogarithm are assumed generic. 

We again start by evaluating the action of the monodromy operator in the left-hand 
side of Eq. (B.4). It is clear that only those terms in the coproduct where the first factor 
depends on a^ contribute a non-zero answer. For generic arguments this implies that the 
polygon associated to this term must contain aj. As an example, for n = 3 and i = 2, the 
term 

"2 




I(a ; ai, o 3 ; a 4 ) <g> I(a 1 ;a 2 ; o 3 ) , 

vanishes under the action of the monodromy operator (.M an+1=a . ® id), whereas the term 

a 2 




ao a 4 

/(oo; 01,02; 04) 8)/(a2;a3;a 4 ) , 

contributes non-trivially to the left-hand side of Eq. (B.4). In Ref. [63] a formula for the 
monodromy of a multiple polylogarithm was given, 

Ma n+1=a J(ao;ai, ■ ■ ■ ,a n ;a n+ i) = 2mI(ao;ai, . . . ,ai_i;oi)/(oi;Oi+i, . . . ,a n ;a n+ i) . 

(B.5) 
Applying this formula to the left-hand side of Eq. (B.4) implies that the monodromy 
operator 'splits' all the polygons that give a non- vanishing contribution at the point Oj. In 
the example n = 3, i = 2 considered above we obtain the splitting 
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ao a4 

[2iril(a ;ai;a 2 ) ■ 1] <8> [1 • I(a 2 ; a 3 ; a 4 )] , 

Summing over the split polygons is equivalent to summing over all pairs of polygons con- 
tributing to the coproduct of the two multiple poly logarithms in the right-hand side of 
Eq. (B.5). Thus, we obtain 

(Ma n+1 = ai <8> id) A(I(a ; ai, . . . , a n ; a n+ i)) 

= (27ri ® 1) A (J(a ; ai, ■ ■ ■ , Oi-iJ a*)) A (/(a*; a i+ i, . . . , a n ; a n +i)) (B.6) 
= A (M an+1 =aJ(ao; ai, . . . , a n ; a n+ i)) , 

which finishes the proof. 
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